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■ ABSTRACT 
I ' 

S ■ We consider the accretion process in a disk with magnetic fields that are 

^ ■ dragged in from the interstellar medium by gravitational collapse. Two diffusive 

processes are at work in the system: (1) "viscous" torques exerted by turbulent 
and magnetic stresses, and (2) "resistive" redistribution of mass with respect to 
the magnetic flux arising from the imperfect conduction of current. In steady 
. state, self-consistency between the two rates of drift requires that a relationship 

Q ■ exists between the coefficients of turbulent viscosity and turbulent resistivity, 

l/^ ■ Ignoring any interactions with a stellar magnetosphere, we solve the steady-state 



O 



equations for a magnetized disk under the gravitational attraction of a mass 
O ■ point and threaded by an amount of magnetic fiux consistent with calculations 

of magnetized gravitational collapse in star formation. Our model mean-field 
^ . equations have an exact analytical solution that corresponds to magnetically 

. diluted Keplerian rotation about the central mass point. The solution yields 

the strength of the magnetic field and the surface density as functions of radial 
position in the disk and their connection with the departure from pure Keplerian 
rotation in representative cases. We compare the predictions of the theory with 
the available observations concerning T Tauri stars, FU Orionis stars, and low- 
and high-mass protostars. Finally, we speculate on the physical causes for high 
and low states of the accretion disks that surround young stellar objects. One 
of the more important results of this study is the physical derivation of analytic 
expressions for the turbulent viscosity and turbulent resistivity. 



Subject headings: stars: pre-main-sequence; planetary systems: protoplanetary 
disks; ISM: magnetic fields; accretion disks; MHD 
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1. INTRODUCTION 

It is universally acknowledged that magnetization is crucial to the accretion mechanism 
in circumstellar disks via the magneto-rotational instability (MRI) first studied in the non- 
linear regime by Hawley & Balbus (1991; see Balbus & Hawley 1998 for a review). Since this 
accretion is the process by which most stars accumulate their masses from the gravitational 
infall of collapsing, rotating, molecular cloud cores (see, e.g., Shu, Adams, & Lizano 1987; 
Evans 1999), and since planets are believed to form from the resulting circumstellar disks 
(e.g, Lissauer 1993; Lin & Papaloizou 1995, 1996; Goldreich, Lithwick, & Sari 2004), a better 
understanding of the mechanism of disk accretion (e.g., Lynden-Bell & Pringle 1974, Pringle 
1981) is desirable for further progress in the fields of star and planet formation. Moreover, 
bipolar outflows and jets are ubiquitous in young stellar objects (YSOs; see Bachiller 1995, 
Reipurth & Bally 2001), while the best contemporary theories for the underlying collimated 
winds intimately involve the combination of rapidly rotating disks and strong magnetic fields, 
either threading the disk itself or belonging to the central host star (Konigl & Pudritz 2000, 
Shu et al. 2000). A major uncertainty in the former models is the strength and geometry of 
the disk magnetic fields. 

In this paper, we consider the global problem of the mass, angular-momentum, and 
magnetic-fiux redistribution in disks around young stars which are threaded by interstellar 
fields dragged in by the process of gravitational collapse and infall. In a future extension 
of this work, we wish to include the interaction of such magnetized accretion disks at their 
inner edges with the stellar magnetosphere generated by dynamos operating in the central 
objects. Such interactions include the loss of angular momentum carried in any outflowing 
wind that develops at the surface of the parts of the disk that rotate sufficiently close to 
Keplerian rotation, a process that is examined in this paper only in terms of whether a wind's 
presence is implied by the prevaihng circumstances. While separate pieces of this problem 
have been attacked by other groups (e.g., Goodson, Bohm, & Winglee 1999; Krasnopolsky 
& Konigl 2002; Kiiker, Henning, & Riidiger 2004; Long, Romanova, & Lovelace 2005), our 
study includes for the first time a perspective that combines analytical calculations with the 
likely levels of magnetic field brought into the disk of a YSO by gravitational collapse and 
infaU (GaUi et al. 2006; Shu et al. 2006). 

Our paper builds on the prescient study of Lubow, Papaloizou, & Pringle (1994, here- 
after LPP) concerning the possibility of disk winds in accretion disks. What distinguishes 
our work from theirs is our concern with the much stronger magnetic fields resulting from 
the process of star formation than assumed by LPP (see §3). Although LPP use a kinematic 
approximation that assumes explicitly an inward drift speed and implicitly a Keplerian ro- 
tation curve, neither assumption turns out to affect the generality of their solutions of the 
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linear induction equation. Nevertheless, when fields are dynamically strong, a solution of 
the full magnetohydrodynamic (MHD) problem is required. Thus, an added component of 
our study is a physical formula for the viscosity associated with the MRI (see §§2.1 and 4.1). 
In other words, while LPP treat magnetic fields as a passive contaminant to an imposed 
accretion fiow, we regard them as underlying the MRI dynamics that drives disk accretion. 
Our own theoretical work on the interaction of an electrically conducting accretion disk with 
a YSO magnetosphere that produces an X-wind and a funnel fiow is incomplete because we 
previously ignored the magnetization of the accretion disk (see the review of Shu et al. 1999 
and the criticisms of Ferreira, Dougados, & Cabrit 2006). 



1.1. Governing Equations 

We wish to calculate the effects of a systematically oriented, poloidal, mean magnetic 
field gathered from the interstellar medium that threads vertically through a circumstellar 
disk that surrounds a newly born star. This field is pinched radially inward by viscous ac- 
cretion through the thin disk driven by the MRI (Fig. 1). To transform the usual 2.5 D 
equations of non-ideal MHD (with steep vertical stratification in z combined with full radial 
dependences in w and axisymmetric motions in the tangential direction Lp) to 1.5 D (inte- 
gration over z) requires that we explicitly treat the vacuum fields above and below the disk 
all the way out to infinity, which is a crucial missing ingredient in all numerical simulations 
of the MRI to date when global fields are present. Fortunately, this transformation can be 
implemented using the Green's function technique used by van Ballegooijen (1989) and LPP 
(see also Shu & Li 1997 and Shu, Li, & Allen 2004, who were unaware of the earlier related 
work on accretion disks until the preparation of the present paper). In 1.5 D, the formulation 
in terms of integro-differential equations is then standard. 

Terquem (2003; see also Fromang, Terquem, & Nelson 2005), made the interesting 
suggestion that toroidal magnetic fields in YSO disks might be strong enough to stop the 
so-called Type I migration of planets and planetary embryos with Earthhke masses. The 
origin of such mean toroidal fields is unclear since they require unclosed z-currents and are 
subject to buoyant vertical loss through the Parker (1966) instability, but similar effects 
could arise for accretion-pinned poloidal distributions. The toroidal fields that arise in this 
paper from the stretching of radial fields by differential rotation vanish in the mean when 
we integrate over z. We shall treat their ffuctuating effects on the turbulent transport of 
angular momentum and matter across field lines as diffusive terms in the non-ideal equations 
of MHD with "anomalous" values for the coefficients of kinematic viscosity v and electrical 
resistivity 77 (see §§2.1 and 4.1). 
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In such a mean-field MHD treatment, the evolution of gas and magnetic field occurs in 
a thin axisymmetric, viscously accreting, disk surrounding a young star that we represent 
as a stationary and gravitating point of mass at the origin of a cylindrical coordinate 
system {zu, (p, z). We denote the disk's surface density by S, the radial velocity of accretion 
in the plane by m, the angular velocity of rotation about the z axis by VL, the component of 
the magnetic field threading vertically through the disk by Bz, and the radial component of 
the magnetic field just above the disk that responds to the radial accretion flow by B^. This 
self-gravitating, magnetized system satisfies the time-dependent equation of continuity, 

as 1 (9 , , , , 

— + -— tuSn =0, 1-1 
ot w ow 



the equation of radial-force balance, 

— w\l = — 

27rS 



GM, + 27r ^ Ko (^) GJ:{r, t) rdr 
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the torque equation, including phenomenologically the effect of turbulent viscous stresses 
(oc z/), 

\ ot ow ) ^ ' w Ow \ Ow ) 

and the induction equation for the vertical component of the magnetic field, including the 
effect of finite resistivity i], 

dBz 1 9 , „ , \ d ( wri ,^ 
—f- + -— (wB^u) = — —Lb^ , (1-4) 

OT W OW W OW \ Zq J 

where, according to Shu & Li (1997), 



oo 



K= I Ko{^)B.{r,tf-^. (1-5) 







wy ' w'^ 



In the vertical averaging over the thickness of the disk to arrive at equation (11-41) . we 
have effectively replaced rj\ocaidB^/dz by its mean value above the mid-plane of the disk, 
TjB^/zo, an operation which defines what we mean by average t]. We shall later consider 
what we mean by the effective half- height of the disk, zq (see Appendix C), but for the 
time being we are content with the intuitive concept. Although it might be mathematically 
more elegant to absorb the combination tj/zq into a single variable denoted, say, by a symbol 
TZ, we retain the more cumbersome notation to keep better contact with the conventional 
microphysics of electrical resistivity. In any case, we assume that zo is much smaller than 
the local disk radius w. 
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In equation fll-2p . the first term on the right-hand side represents the mass per unit 
area, S, divided into the radial component of Lorentz force per unit area due to magnetic 
tension, J^i?^/c, where is the current density integrated over the thickness of the disk 
and equals ci?+/27r by Ampere's law. The second and third terms on the right-hand side 
represent the contributions to the force per unit mass associated with the stellar gravity of 
point mass M^, and the self-gravity of the gas of surface density S in the disk. To lowest 
order in the aspect ratio, zq/w <^ 1, we have neglected the pressure forces of the matter 
and the magnetic field (see Shu & Li 1997). The important astrophysical point is that the 
centripetal acceleration on the left-hand side of equation (11-21) is not given a priori but arises 
in response to the forces of (1) magnetic tension, (2) stellar gravity, (3) self-gravity of the 
disk gas, and (4) gas and magnetic pressure forces. We have ignored (4) explicitly, and 
we shall presently ignore (3) also. The rationale is that (3) and (4) are generally small in 
comparison with the stellar gravity term. Their inclusion would only yield small corrections 
at the expense of rendering the resulting problem intractable except by numerical attack. 
In contrast, the magnetic tension force is always present at a level dictated by the amount 
(and distribution) of magnetic flux threading the disk. 

The kernel in equation fll-51) is given by 

The function -ft'o(0 is plotted in Figure 1 of Shu, Li, & Allen (2004) with mathematical 
properties described in their Appendix A. It has the asymptotic behaviors: 

iro(O) = 1 and Ko{i) --^ as ^ ^ +oo. (1-7) 

2^^ 

On a microscopic level, field diffusion in lightly ionized gases involves, in principle, 
three non-ideal effects (see Wardle & Ng 1999): Ohmic resistivity (because ions are knocked 
off field lines by collisions), the Hall effect (because electrons move differently than ions 
under electromagnetic fields and collisions), and ambipolar diffusion (because neutrals do not 
feel electromagnetic forces directly, but are subject to them indirectly because of collisions 
with the ions). For an axisymmetric problem, the Hall term vanishes, and the remaining 
two effects can be accorded the same treatment by defining an effective resistivity r/ given 
by the sum of the Ohmic resistivity ?7ohm and the contribution from ambipolar diffusion 
tB'^/Attp, where p is a representative local volume density and r is the mean collision time for 
momentum exchange between a neutral particle and a sea of charged particles or particulates 
(if charged dust grains are important): 

'7 = ^Ohm + X^. (1-8) 
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For mean-field magnetohydrodynamics, we assume that the same relation applies phenomeno- 
logically, and henceforth, when we speak of the "resistivity," we mean the effective value fll-Sp . 
appropriately generalized to include turbulent fluctuations. In Appendix A, we comment on 
some alternative formulations of the turbulent diffusive processes that result in the same 
conclusions when the system is in steady state. 



1.2. Steady State 

In steady state with a spatially constant mass-accretion rate Md, the equations of con- 
tinuity, torque, and advection-diffusion of magnetic field simplify to 

zuTjU = const = -, (1-9) 

27r 

Mdzu^n = -27rro^Si/^ (1-10) 

aw 

BzU = -B,- = , 1-11 

ZTTWh Zq 

whereas the centrifugal balance and the radial magnetic field just above the disk plane 
are still given by equations fll-2p and fll-51) with no modifications except that there is no 
time dependence in the arguments of S and Bz- Elimination of from equations (11-101) 
and (11-111) then yields the self-consistency requirement, 

B+ u dn , , 

rj ^ = . (1-12) 

zqBz dzu 

In choosing the integration constants as above, we are implicitly allowing the origin (the 
star) to be a sink for mass but not for magnetic fiux (or for angular momentum). To make 
this clearer, multiply equation (11-41) by 27rct7 dw and integrate in radius from the origin to a 
position just outside the star i?+ (which we will ultimately let — > 0"*"). The result yields 

-2rrK(BzU+^-^) , (1-13) 

x7 = _R+ 



dt \ Zq 



where is the magnetic fiux accreted by the star. The assumption that equation (11-111) 
holds outside the star, in which the radial advection of magnetic field is everywhere balanced 
by the radial diffusion associated with diffusive effects, then implies that 



(1-14) 
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i.e., all the magnetic flux that is brought in by the gravitational collapse involved in star 
formation is contained in the disk, if the star plus disk forms a closed system. We emphasize 
that equation fll-14p is an assumption, not a deduction from first principles. It is astronomical 
observations, not theory, that tell us that little interstellar flux is brought inside stars. Indeed, 
Mestel & Spitzer (1956) already recognized more than fifty years ago that stellar fields would 
measure in the megagauss range if even a small fraction of the original interstellar flux were 
to appear on a young star's surface (see also Shu et al. 2004, 2006). 

In this context, it is important to realize that magnetic reconnection will only destroy 
field lines that close within the disk. Indeed, the annihilation of field loops is essential for 
fluid transport through the disk, as discussed in §4.1 below. In contrast, the open, mean, 
magnetic fields pictured in Figure 1 have roots that extend to the interstellar medium. These 
field lines are responsible for the flux (integral of Bz over the area of the disk), and this flux 
cannot be lowered without modifying the interstellar currents (assumed to vanish at all 
finite distances in the derivation of eq. [l-5j ). for example, by arbitrarily adding spatially 
large loops of current that yield a flux of opposite sign to the original interstellar valued 

In what follows, then, we choose to write equation (11-121) as 

5+ = a^Bz, (1-15) 

where we define the dimensionless auxiliary variable a by 

^ _ zqu f w dn\ 

a = — — . 1-16 

zur] \\l am J 

For realistic configurations, we expect ~ B^, i.e., we anticipate that a is an order unity 
quantity. Since wQ~^dQ/dw is also of order unity, the resistivity rj must be smaller than the 



^One might still worry about the fate of open field lines with a fixed flux that connects a rapidly rotating 
disk with a slowly rotating interstellar cloud. Wouldn't such field lines get twisted up in time and torque 
down the disk? The answer is yes, but for disturbances traveling a poloidal distance As along a field line out 
of the plane of the disk, the change in azimuthal angle Aip experienced by a loaded field line is approximately 
given by the disk rotation rate times the Alfven crossing time across the region: A(p ^ ^.{As/va)- During 
the infall stage, the mass loading is so considerable that the ability even to form a disk is in jeopardy (see, 
e.g., Galli et al. 2006). After infall has stopped, the Alfven speed is very large near but above the disk, 
and it decreases as the field line begins to penetrate the interstellar cloud. In this stage, all the twist is 
at the cloud end; virtually none is at the disk end. In other words, equation p-25p below will hold to 
good approximation. There is little continued torquing of the disk because the rate of angular momentum 
transport per unit area by torsional Alfven waves is limited by the angular-momentum density pzu'^fl times 
the Alfven speed va = B/y/Anp, with the small p above the disk of the former more than canceling the 
in the denominator of the latter. 
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(turbulent) viscosity u hj a. factor given roughly by the local aspect ratio of the disk, i.e., 

^ ~ ^ « 1. (1-17) 

Equation fll-17p is in agreement with the assertion in equation (39) of LPP that the di- 
mensionless ratio of interest is not rj/v but [t] / v){w / Zq) , which must be of order unity 
for magnetic fields to be bent by an order unity amount from the vertical direction. The 
contrary assertion by Riidiger & Shalybkov (2002) arises because they arbitrarily assume 
a uniformly-rotating halo of conducting matter in which the differentially rotating disk is 
embedded. Some aspects of the problem that Riidiger & Shalybkov consider could apply 
to the interaction of the inner edge of an accretion disk with a uniformly rotating stellar 
magnetosphere (e.g., Kiiker et al. 2004). The latter is the context of X-wind theory. In 
the next subsection, we give an explicit justification of the order-of- magnitude estimate in 
equation fll-17p when the disk mass is small enough to allow us to ignore its self-gravity. 



1.3. Exact solution when self-gravity of the disk is neglected 

If we may ignore the last term (disk self-gravity) in equation (ll-2p in comparison to the 
other terms on the right-hand side (magnetic tension and stellar gravity) , then the condition 
of centrifugal force balance becomes 

where is given by equation (11-50 . Consider now the important case when f2 is appropriate 
for a thin disk in quasi-Keplerian rotation. 



^ = f\-^] ' (1-19) 



where / is a constant less than 1 because of partial support of the disk against the stel- 
lar gravity by the magnetic tension of the poloidal magnetic fields that thread through it, 
brought into the system by the process of star formation. In order for these conditions to 
be mutually compatible, equation fll-18p requires 

B, = [27r(l - f)GM,w-^J:] . (1-20) 

The substitution of equation fll-20p into equations (ll-15p and (11-50 now results in a hnear 
integral equation for S^/^ when a is known, 

aivj)J:'/^{w) = r kJ -] a-\r)i:^l\r) —. (1-21) 

In \W / ZU 
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Because S^/^ enters linearly into equation (11-211) and ot enters nonlinearly, the above 
relationship yields a constraint on the allowable solutions for both S and a. Consistent with 
equation fll-19p . suppose, for example, that S is given by a power law: 

= Cw~'^\ (1-22) 

where C and £ are constants. Then, equation (11-211) requires to be a positive dimensionless 
constant: 

/■oo 

o? = h= \ KoiOr'd^- (1-23) 







The integral has a finite value for i between —2 and 1 with Iq = 1. Table 1 gives a 
tabulation of numerical values for the astronomically interesting range of i from to 1 
where S is a declining function of radius. The last row shows the inclination angle i that the 
surface field makes with respect to the vertical direction as computed from tani = B^/B^ 
and equation (I1-25P below. 

Equations (I1-23P and (11-190 now allow us to deduce from equation (I1-16P the required 
relationship between the resistivity and viscosity as 

" - ' ^^^) . (1-24) 



p 21 I \vj J 

which supplies the missing numerical coefficient to our previous estimate (I1-17P . The Prandtl 
combination rj/u is required in a steady-state accretion disk to have the specific ratio (11-241) 
if u = —?)v/2w arising from the viscous transport of angular momentum (cf. eqs. |l-9] and 
|l-10j ) is the same drift speed needed for the resistive diffusion of matter across stationary 
field lines, u = —{ri/zo)B^/Bz (cf. eq. |l-llj ). The two formulae for the drift velocity 
express succinctly why rj only needs to be a small fraction of u: unlike viscosity, resistivity 
is not acting to mix quantities on a large scale of w, instead, it is trying to annihilate the 
oppositely-directed mean radial-fields 5^ on either side of the mid-plane distributed on a 
small scale zq. If we consider the diffusivity associated with the radial distribution of the 
current (Appendix A), then that diffusivity rjj is approximately equal to u. 

The corresponding relationship between the radial component of the magnetic field at 
the upper surface of the disk and the vertical component at the mid-plane is given by equation 
(f^T5|) as 

Bt = hB.. (1-25) 

The vertical field is itself given by equation (ll-20p when we know the surface density from 
equation (ll-22p . i.e., 

B, = J7'/' [27r(l - f)GM,C] w'^^+'l (1-26) 



Table 1. Values of le and inclination angle i 



e 0.1 0.2 1/4 0.3 5/16 3/8 0.4 1/2 0.6 0.7 0.8 0.9 1 
le 1 1.149 1.326 1.428 1.542 1.573 1.742 1.818 2.188 2.726 3.598 5.304 10.34 oo 
i 45° 49.0° 53.0° 55.0° 57.0° 57.5° 60.1° 61.2° 65.4° 69.9° 74.5° 79.3° 84.5° 90° 
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To obtain the identification of tlie coefficient C, we note tliat equation (ll-lOp implies 

Sttu 



1-27) 



Tlius, tlie adoption of fll-22p is an implicit assumption that the viscosity varies as a power-law 
of zu in steady state given by 



.^1^1^- fl-28) 

To make further progress, we need to have a physical theory for the kinematic viscosity 
u (see §§2.1 and 4.1). This will allow the last unknown quantities i and C to be eliminated 
from our power-law solution for Q, S, and B^. We discuss first, however, in §§1.4 and 1.5 
two results that hold for more general diffusivities. 



1.4. Disk winds 

There are two criteria necessary for a cold wind to be magnetocentrifugally driven from 
the surface of a rotating disk (cf. the review of Konigl & Pudritz 2000). Without detailed 
specifications of the physics of the viscosity or the resistivity, equation (11-251) allows us to 
confirm the conclusion reached by LPP concerning the first criterion: 

so that the footpoint field would then bend with an inclination angle i from the vertical by 
more than 30° (Chan & Henriksen 1980, Blandford & Payne 1982). Table 1 shows that the 
above criterion is comfortably satisfied for any disk where the surface density declines with 
radius, i > 0. In particular, i = 1/4 in equation (I1-26P corresponds to the famous case 
Bz oc zu~^^^ considered by Blandford & Payne (1982) and yields i = 55.0°. LPP state their 
wind-launching criterion in the form that [r]/ v){w / zq), which equals 3/(2/^) according to 
equation (11-241) . should be less than 1.52a/3. This result is almost identical to the criterion 
(ll-29p . LPP's calculation is for a finite disk with nonzero inner and outer radii, embedded 
in a background field of uniform strength pointing in the vertical direction, whereas our 
calculation is formally for an infinite disk with a trapped interstellar flux. The negligible 
astrophysical difference between 1.52 and 3/2 implies that none of these idealizations matter 
to the first criterion for driving a disk wind. 

Unfortunately, the satisfaction of the magnetic criterion (11-290 by itself is not a sufficient 
condition for the appearance of a significant disk wind. Equation (12-121) . derived from the 
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consideration of vertical hydrostatic equilibrium, shows that if the fractional departure from 
Keplerian rotation 1 — is large in comparison with the aspect ratio, A = zq/w divided by 
then the square of the characteristic thermal speed in the disk interior, a?, is given by 

a^-|A(l-f)^. (1-30) 

On the other hand, in order to drive a disk wind, the square of the thermal speed at the disk 
surface, a^, must be greater than a fraction (say, 1/4) times the virial imbalance between 
the gravitational potential and twice the specific kinetic energy in disk rotation: 

al>\{l-f)^. (1-31) 

The 1/4 on the right-hand side has the following approximate justification. Parker's solution 
for a thermally-driven spherical wind gives one factor of 1/2 (see Parker 1963); this 1/2 
becomes 1/4 because a particle in Keplerian rotation already has 1/2 of the energy needed 
to escape. No other factors are then included because magnetocentrifugal effects do not help 
MHD winds in making the sonic transition (see the discussion of Shu et al. 1994). 

Except for the effects of heating by external irradiation, the square of the thermal speed 
al of the gas at the disk surface is likely to be small in comparison with its value in the disk 
interior when the disk is vigorously accreting, the condition needed to allow equation (12-121) 
to hold. The inequality (11-311) is then inconsistent with equation (ll-30p . implying that strong 
disk winds cannot be driven when 1 — P much exceeds the small number A (the disk aspect 
ratio zq/w). Otherwise, the surface that corresponds to smooth slow-MHD crossing would 
lie at so many scale heights above some nominal disk surface that the associated mass-loss 
rate would become negligibly small. LPP avoided this problem by their implicit assumption 
that the fields threading the disk were weak and therefore had no effect on the disk's assumed 
Keplerian rotation. Wardle & Konigl (1993) examined the same issue in a /oca/ treatment of 
the launch region for disk winds assuming ambipolar diffusion to be the physical mechanism 
that loads field lines. They reported that the effect is present, but compensating factors 
exist that allow wind mass-loss rates to be a small fraction of the disk accretion rate, with 
a very sensitive dependence on the ratio of the orbit time to the ion-neutral coUisional time 
r (see Fig. 12 of their paper). The global treatment given in this paper, which includes an 
assessment of the likely levels of magnetic field strength to result from the process of star 
formation, indicates that the problem is more severe, perhaps even insurmountable for FU 
Orionis and T Tauri stars, although the situation may yet be rescued for the outer disk 
regions of embedded low- and high- mass protostars where A is not so small (see §3). 

Font et al. (2004) propose that photo-evaporation is a more likely source of the slow, 
warm disk-wind observed by Kwan et al. (2006) in T Tauri stars. A photo-evaporative wind 
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could reach a higher terminal velocity, but not a larger mass-loss rate (limited by the X-ray, 
EUV, or FUV photon-flux reaching the surface of the disk at radii of a few AU or greater), 
because of the boost to the gas by magnetocentrifugal fling after the sonic transition is made. 
The combined effect could lead to a signiflcant loss of angular momentum from the system 
in the late stages of YSO evolution that is not taken into account in our treatment. In other 
words, outer disk winds may realistically arise even if 1 — is not small. 



1.5. Viscous and resistive dissipation rates 

In steady state, the net emergent radiation from the upper and lower surfaces of a thin 
disk, after accounting for the irradiation of the central star, has to carry away the sum of 
the energies generated by viscous and resistive dissipation, whose rates per unit area are: 

{--B.) ^ (1 - f )^^ 

2tt J \ c / Znw^ 



Y^J,E,= {-^ { —B, = (1 - r)^-^, (1-33) 



In the above, = ci?+/27r and E^p = —uB^/c are, respectively, the mean current density 
(integrated over z) and electric fleld in the ip direction in the rest frame of the plasma tied 
to mean Bz relative to which the bulk of the matter in the disk is drifting at radial velocity 
u = —Md/27izuT, = —?)u/2w = —{rj/zQ)B^/Bz. To perform the last step in equation f ll-33p . 
we have used equations fll-lSp and (11-191) to eliminate B'^Bz/2-n and VL. 

An alternative expression, Y = ri{B^Y /2t:zo, makes more apparent that Y represents 
the resistive dissipation, which feeds on the magnetic tension. In contrast, the viscous 
dissipation \l/ feeds on the disk shear. In the former case, a Lorentz force drives electric 
currents that generate heat by friction between the various charged and non-charged species; 
in the latter, heat is generated by fluid elements "rubbing" tangentially against each another. 
We speculate that a fraction of the energy released by resistive dissipation in the disk may 
go into accelerating suprathermal particles that give proto-planetary disks higher ionization 
rates than conventionally estimated (cf. Goldreich & Lynden-Bell 1969 vs. Sano et al. 2000). 

The coefficient 3/^/2 in equation (11-321) differs from the standard result by the factor p 
because the rotation law is only quasi-Keplerian, and implies a local rate of energy release 3 
times as great as one might have expected from the loss of orbital energy because of accretion 
(see, e.g., Lynden-Bell & Pringle 1974). The difference is made up by a viscous torque that 
transfers energy from the inner disk to the outer disk, a debt that has to be repaid if we were 
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to examine the details of the interaction of the disk near its inner edge with a star of finite 
size and, perhaps, magnetization. For example, if one applies a zero torque condition at an 
inner boundary corresponding to a stellar magnetopause that corotates at the angular 
rate n(i?x) = f{GM^/ R^^/'^ , then standard arguments (Lynden-Bell & Pringle 1974) show 
that the total viscous energy dissipated by bringing matter through our accretion disk from 
infinity to i?x equals 

i/'^. (1-34) 

In contrast, if we multiply equation fll-33p by 27r zudzu and integrate from R^ to oo, we 
derive that the rate of resistive dissipation of energy in the disk equals 

(l-f)^. (1-35) 

The sum of the viscous and resistive dissipation rates, equations (11-341) and fll-35p . equals 

1 GM.Md 



1 - ^f'j (1-36) 

which is not the total rate of gravitational potential energy release, GM^Md/R^, because 
an amount R\VL^{R^/2 = pGM^,/2Rx is still retained by each gram of disk matter at the 
disk/stellar- magnetopause boundary as specific orbital energy. In actual practice, as will 
be shown in a future publication, / is raised to unity as the latter boundary is crossed by 
the swing of B^Bz from positive to negative values and inner-edge effects (assuming the 
magnetopause is not squashed by the accretion flow to the stellar surface in quasi-steady 
state), so GM^/2Rx of specific energy in the accreting matter is available for budgeting in 
a funnel flow or X-wind (Shu 1995). Notice that in this description, electromagnetic fields, 
although responsible for the microphysics of viscous and resistive dissipation, act on the 
macroscale merely as catalysts for converting gravitational energy into other forms. These 
other forms, in steady state, involve no change of the magnetic energy because the magnetic 
fields have been assumed to remain constant in time. 

Let us compare the expressions (ll-32p and (11-330 which hold at radii w ^ R^. Then 
we easily calculate that heat generation by viscous dissipation dominates over resistive dis- 
sipation when / > a/2/5 = 0.6325. The resistive contribution is typically not negligible; for 
example, it is 37.5% of the viscous contribution when / = 0.8. In spirit, if not in detail, 
our ideas then follow those of Lynden-Bell (1969), and we can imagine "resistive accretion 
disks" as well as their "viscous" counterparts (cf. the FU Orionis model of §3). 
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2. VISCOSITY ASSOCIATED WITH MRI 

Returning to our quest for the viscosity v to be used in our power-law solution in §1.3, 
we would like to benefit from the many numerical simulations that have been performed 
of the MRI since the pioneering studies of Balbus and Hawley. However, few experiments 
have been done that are of direct relevance to the problem in star formation that we address 
in this paper. Most simulations miss one or the other of the crucial ingredients of being 
both global and having nonzero net flux. When field lines extend to infinity (necessary to 
have nonzero net flux), rather than close within the system, consideration of the behavior 
of the fleld within a few vertical scale heights of a spatially thin disk suffices only if one has 
included knowledge of what those fields do at large distances. The application of boundary 
conditions at smaller distances will generally exert extraneous stresses. Because of these 
difficulties, many simulations are both local and have zero net fiux, in that a small portion 
of a shearing sheet or layer is threaded by mean vertical magnetic fields inputted initially 
to vary sinusoidally in the radial direction. Radial mixing and reconnection can destroy 
most of the initial vertical field in such simulations, so that the turbulent state reached 
asymptotically in time is largely independent of the assumptions of the "initial" state. The 
finely resolved study by Silver & Balbus (2006) does include the effect of a systematically 
directed field of a single sign, but their simulation is not global (and therefore does not 
develop a B^ comparable to B^), and it is performed for a gas pressure 800 times larger 
than the magnetic pressure, i.e., with implied magnetic fields that are too weak to be useful 
for our studies here. 

Stone et al. (1996) and Miller & Stone (2000) performed, to our knowledge, the only 
well-known simulations of a thin disk threaded by a systematic, large-scale, nonzero, vertical 
field Bz- The computations were semi-global in spanning a larger than usual, but still 
hmited, range of z. (The work by Fleming, Stone, & Hawley 2000, assumes periodicity 
in the z direction, which does not faithfully represent the dynamics of a thin disk.) The 
thin-disk cases with non-zero net fiux behave completely differently from the other more 
frequently studied configurations, whose initial states have only toroidal fields, or, at least, 
zero average Bz- The systems in the simulations of Stone et al. (1996) and Miller & Stone 
(2000) quickly become magnetically dominated, unlike the usually considered circumstance 
where the gas pressure is much greater than the magnetic pressure. The rapid evolution then 
prevented the authors from examining the astrophysical consequences of the configuration 
most likely to be relevant to investigations in star formation. 
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2.1. Turbulent Viscosity 

In the absence of relevant numerical simulations, we give the following order of mag- 
nitude argument on the basis of mixing- length ideas (Prandtl 1925). When a radial field 

is present in a field of differential rotation, we expect that field to be sheared and yield 
an azimuthal component B^^. The tendency for electrically conducting fluids to flow along 
the field direction suggests that fluctuations in the radial velocity 6u will be related to the 
horizontal fields and the shear rate by 

B^5u ~ B^zu- — 5zu, (2-1) 

where Svu is the radial mixing length and has the same sign as 6u. Notice that the induced 
Bi^ has systematically the opposite sign as B^ if dVL/dw is negative. The systematics of B^ 
relative to B^ lead to the desired "viscous" torque. 

Differentially rotating fluid parcels displaced from their equilibrium positions that pre- 
serve their specific angular momentum gyrate in epicycles about a guiding center character- 
ized by an epicyclic frequency k that numerically equals in a quasi-Keplerian disk (Binney 
& Tremaine 1987). Although other forces are also at play in a magnetized accretion disk, we 
assume that mixing-length scales of greatest interest for the transport of angular momentum 
have a correlation time between 5w and 5u that is similarly given by ~ i.e., 

5u ~ Vl5w. (2-2) 

Equation 02- ip can now be written 

„ „ vo dVt 

B^^B^- — . (2-3) 
\l dzu 

The component of Maxwell stress responsible for exerting torque, B^Bi^/Att, integrated over 
the thickness of the disk, can then be approximated as 

(B+fwdVL , , 

27r \l dw 

where JF is a form factor that comes from the vertical integration, and that also corrects 
for all the order-of- magnitude approximations used to arrive at this point. If the term (12-41) 
is the dominant contribution to the "viscous" stress modeled in equation (11-31) . then the 
associated "kinematic viscosity" equals 
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A more picturesque "derivation" of equation fl2-5p which explains why the correlations 
do not involve quadratic products of fluctuating quantities (they actually do) and why the 
mixing length 5w seemingly dropped out of the calculation (it should not) is given in §4.1. 
We proceed here to build confidence in the case by first demonstrating that the adoption of 
equation (12-51) leads to reasonable astrophysical results. 

To follow the flux redistribution, our problem is formulated using rather than 
so in steady-state we set S+ = UB^ (cf. eq. |l-25] ) and get 

where D is a dimensionless coefficient given hy D = IjT. Although T is the more fun- 
damental quantity, as we shall see in a later study of the interaction of disks with stellar 
magnetospheres, we shall use D in this paper as the relevant dimensionless parameter to 
obtain v from observations and simulations (see the discussion of Appendix B). As long as 
T is not too small compared to unity, D is an order unity quantity, provided the entire disk 
layer undergoes vigorous mixing from the MRI. 

With V given by equation (l2-6p . equation (I1-27P yields the following expression for the 
vertical magnetic field: 



where we have used equation fll-19p to express the angular rotation rate VL and where A = 
zo/zu is the aspect ratio of the local disk height to disk radius. From equations (11-200 and 
(I2-I2p . we are then able to recover the surface density as 

f f h \ Ma . s 



1-/2 \?,TxDAJ {GM.wy/^' 

Equation (l2-6p holds with nonzero D only as long as (1) good magnetic coupling exists, 
and (2) the criterion for the MRI instability is satisfied, that the magnetic pressure be smaller 
than the gas pressure. If we define a fiducial square of the thermal velocity a'^{w) by the 
gas pressure at the midplane, P{w,0), divided by the characteristic volume-density in the 
disk, S(cc7)/2zo(^), equations (12-71) and (12-81) imply that the ratio of the magnetic pressure 
to gas pressure at the disk's midplane is then given by 

B!zo ^ (1 - P) ( AG MA 

47rEa2 2/, \ a^zu J ' ^ ' 

where we have expressed zq = Aw. On the other hand, analysis of the vertical hydrostatic 
equilibrium of the disk using the method of Wardle & Konigl (1993; see also Ogilvie 1997 
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and Shu & Li 1997) gives (see Appendix C): 

1 r. „o. . .01 CM, 



[/,(1-/2)A + A2]^. (2-10) 



In principle, a proper physical treatment would require us to obtain c?{w^ by computing the 
volumetric viscous and resistive heating as a function of z and balance it against the heat 
transported by radiative transfer and thermal convection vertically out of the disk. Once we 
have obtained a^{w\ we could then solve equation (12-101) as a quadratic equation for the 
disk aspect ratio A. Such a physically involved treatment is beyond the scope of the present 
paper, and for the astronomical and pedagogical sake of obtaining numerical examples, we 
shall assume the luxury of specifying Aiw^ semi-empirically as a power-law (see §2.3). 

The justification for approximating A as a power-law follows. If disks are spatially thin, 
A is small compared to unity. There are then two different regimes of physical interest. The 
first case arises when the departure (1 — p) from Keplerian rotation is small, and the first- 
term on the right-hand side of equation (12-lOp is negligible in comparison with the second 
term, resulting in the approximation: 



In case (12-111) . the contribution of the magnetic pressure is ignorable for the vertical hy- 
drostatic equilibrium, and the draw of the stellar gravity toward the midplane keeps a cool 
accretion disk spatially thin. The second case arises when the departure from Keplerian 
rotation (1 — /^) is not small, and the first-term on the right-hand side of equation (12-101) 
dominates over the second: 



, 2 
A^ — 
If 



a?'UJ 



for 1 - /2 > ^. (2-12) 



In case (12-121) . the disk is kept spatially thin, not by stellar gravity, but by the inward press 
of the component of magnetic pressure B^/Stt which increases outward from the midplane 
z = 0, where its value is 0, to the surface, where its value is {B^^Y/Sn = IfB'^/SiT. 

From the community experience in disk thermal-physics, it is well-known that numer- 
ical solutions frequently show two types of power-law solutions for the vertically-averaged 
temperature oc a^{zu), namely, oc for passive, irradiated disks, and oc w^^^'^ for 

active, accretion-powered disks. For oc tu"^/^, case (12-1 ip yields A(zu) oc cc^/^ and case 
(12^ yields A{w) oc w^/^. For oc ro^^/^ case f lCT]) yields A{w) oc zu^/^ and case (|2^ 
yields A{w) oc ro^/^. Therefore, a power-law description for disk flaring, A^tu) oc tu", with 
n = 1/4 typically or n = 1/8 or 1/2 at the extremes, has a theoretical basis. 



Apart from the issue of sufficient ionization, the requirement for MRI being present in 
the mid-plane is that the left-hand side of equation (12-91) should be less than unity. This 
requirement is automatically satisfied for all our disks because the substitution of equation 
fl2^ into the right-hand side of equation ([23]) gives a value (1 - - + Ah], 

which is always smaller than 1 for Ii > 1, i.e., for i > (Table 1). In particular, in the 
regime where equation (I2-I2p holds, u from equation (12-61) becomes 



which has the same form as Shakura & Sunyaev (1973) viscosity, u = assd"^/^ with ass = 
2D /If. The strongly magnetized disks of this paper therefore both automatically satisfy the 
criterion that MRI exists in the midplane and have equivalent Shakura-Sunyaev alpha's of 
order unity if D ~ 1 . By compressing the midplane gas density and pressure to higher values 
than gravity can achieve alone, such disks always operate at nearly the maximum efficiency 
for viscous transport, if D ~ 1, without shutting down the MRI. Aficionados of MRI like to 
say that it is present in thin disks for arbitrarily low levels of magnetic field; now they can 
add that it is present for arbitrarily high values too, provided Ii = B^/Bz > 1. 

For later reference in discussions of disk fragmentation, we record that the local dimen- 
sionless mass-to-fiux ratio in the disk is given by (see Nakano & Nakamura 1978, Basu & 
Mouschovias 1994, Shu & Li 1997, Krasnopolsky & Gammie 2005): 



The supercritical condition A > 1 is necessary, but not sufficient, for local disk fragmentation. 
We must also examine, at least, the Toomre (1964) Q parameter, which for gaseous disks 
must be less than unity for local gravitational instability. Thus, also for later reference, we 
note that associated with equation (12-101) is a Toomre Q, which is given by the formula 
Q = Qa/TrGH for a disk in quasi-Keplerian rotation. In principle, a should be computed 
from the considerations of heat balance outlined earlier, but for fiducial purposes, we use the 
value of a associated with case (I2-I2p : 



In what follows, we take the combination, A > 1 and Q < 1, as necessary indicators for local 
gravitational instability. With a "standard" fiaring law, A oc w^^'^, both criteria favor the 
outer regions of a disk for the possible occurrence of disk fragmentation. 




(2-13) 




(2-14) 




(2-15) 
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Finally, for arbitrary values of /, equations (11-271) and fll-24p give the viscosity and 
resistivity generated by the MRI instability in the disk as 



AD{1 - f 



2\ 



{GM,wY'\ (2-16) 



fh 

ri = ^-^-^^jf^iGMM'^'- (2-17) 



Apart from the factors involving A, D, and /, these expressions show that the natural 
scaling for both diffusivities is the specific angular momentum of the matter in Kepler orbits, 
{GM^^zuY^'^. Notice especially the lack of any parametric dependence on the assumed mass 
accretion rate M^. 



2.2. Enclosed disk mass, magnetic flux, and angular momentum 



To make equations (12-71) and (I2-8P consistent with equations (I1-26P and (11-221) . the 
product DA has to be a power law of zo. We adopt the natural assumptions that D equals 
a constant and the disk flares as a power-law, so that the aspect ratio A{r) at radius r is 
related to its value A{w) at radius w by the formula: 



A{r) = A{zu){r/w) 



(2-18) 



where n is the flaring exponent, which is positive deflnite if shadowing does not occur. Had 
we adopted these assumptions from the start, together with the hypothesis (l2-6p . we could 
have shown that not only are the discovered power-law solutions possible in steady state, but 
they are unique. The relationship between the exponent n and the exponent i in equation 
(I1-22P can be found from equation (I2-8P with A{zu) oc zu"', namely, —2i = —n — 1/2, or 



;i + 2n)/4. 



(2-19) 



For later reference we note from Table 1 that Ii = 1.573, 1.742, or 2.188 for disks with low, 
typical, or high power-law flaring, n = 1/8 and i = 5/16 (i.e., S oc ti7" 

-5/8)^ n = 1/4 and 

i = 3/8 (i.e., S oc zu^^^^), or n = 1/2 and £ = 1/2 (i.e., S oc zu^^), respectively. 



We wish now to compute the enclosed mass in the disk inside a radius zu, 



Jo 



Alp 



_3(3-2n)J VI-/ 



/ 



_DA{w) 



(2-20) 



Thus, M(i{w) is a multiple of the mass that accretes through the disk during a Keplerian 
rotation period at that radius. This multiple depends on the combination {DA)^^ and on 
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the effective value of /. Except for a numerical factor of order unity, we easily verify the 
interpretation that the enclosed mass Ma{w) results from accretion at a rate Md over a 
viscous diffusion time scale w'^/u. 

According to the discussion in §1.2, the magnetic field brought in by infall is contained 
as magnetic flux threading the disk. Inside any radius w where the disk is in steady state 
the enclosed magnetic flux is given by 



$d(ro)=27r/ B,{r)rdr = 27c 
Jo 



4 



3 - 2n 



2/ 



3DA{va) 



1/2 



iGM^M^wY^^ (2-21) 



If the disk's mass is negligible in comparison with the star's, the system's dimensionless 
mass-to-fiux ratio enclosed inside vj equals 



27rGi/2M, 



2n 



3DAiwy\ / GMl 
2/ J \M~ 



1/4 



(2-22) 



Infall models with field freezing until small radii yield A=k ~ 1 to 4 (Galli et al. 2006), 
consistent with the polarization findings of Girart et al. (2006) of an hourglass shape in NGC 
1333 IRAS 4A. Field slippage during the collapse reduces the enclosed flux for a low-mass 
protostar by a further factor of 2 to 3 at the radii ~ 300 AU that their disks are likely to 
occupy (see Fig. 3 of Shu et al. [2006] when -Rohm in that paper has a value ~ 10 AU). Thus, 
in §3, we shall adopt an enclosed mass-to-fiux value of Aq = 4 for the system as a typical 
outcome of the star and disk formation process. If insufficient time has elapsed for the disk 
to reach steady state inside the radius where A* = Aq, investigations of the affected regions 
should make use of the time- dependent equations with which we began this paper (§1.1). 

It is extremely informative to compute the enclosed mass at a radius -R$ where A*(i?$) = 
Ao, i.e., at a radius where the disk contains the entire flux brought in by star formation: 



For the typical case, n 



1/4, £ = 3/8, and h 



(3 - 2n)h 



l_/2- 

1.742, equation fl2-23p becomes 



(2-23) 



0.5444\ M, 



^0 



(2-24) 



For a closed system in which infall has ceased, so that Aq remains a fixed constant, disk 
accretion must decrease Mrf(i?$) relative to M*, and therefore the departure from Keplerian 
rotation, (1 — /^), must grow with time. This trend arises because viscosity drains mass 
from the disk onto the star, while resistivity can only cause the redistibution of flux within 
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the disk but cannot change the total, making the specific magnetization (inverse A) rise 
with time. In §3 we shall combine equation (12-241) with an assumption of a disk's "age" to 
estimate the numerical value of the important parameter /. 

To appreciate concisely the dynamical consequences of disk magnetization, we note that 
equations (12-141) and (12-221) imply the interesting reciprocity relationship: 

A(-)A.(-) = (^) (^) , (2-25) 

where the right-hand side is a constant that depends only on / and C. = {1 + 2n)/4. Except 
for such constants, a similar relationship was derived by Shu, Li, & Allen (2004) in their 
analysis of magnetic levitation of pseudo-disks by strongly magnetized central objects (see 
their eq. [65]). Although disks differ from pseudo-disks in being (partially) centrifugally 
supported, and although the inner parts of a magnetized disk differ from a split monopole in 
their detailed interaction with the outer parts of the same magnetized disk, the principles are 
qualitatively similar and provide physical insight into why such magnetic support (levitation 
of the outside by the inside) causes the rotation to occur at sub-Keplerian rates. 

Finally, if the mass is mostly in M^,, the enclosed angular momentum of the parts of the 
disk that are in steady-state accretion is given by 

Jd = 27, nr)r^n{r)rdr = f (j^^ Md{w){GM,wf\ (2-26) 

Equations (12-211) . (12-201) . and (12-261) show that the enclosed disk flux, mass, and angular 
momentum represent a sequence of decreasing central concentration. 



2.3. Semi-numerical formulae 



For the convenience of the reader, we express the results of the analytical theory in the 
following semi- numerical form (assuming n = 1/4): 



B,{w) = 8.89 X 10-3^-1/2 f^^^ 
^ ^ VO-5 Mq/ 



1/4 



1/2 



X/ 



1/2 



0.1 

A{w) 
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G, 



(2-27) 
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(2-32) 
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(2-34) 



(2-35) 



3. ASTRONOMICAL EXAMPLES FROM STAR FORMATION 



To give astronomical context to the theory developed so far, we consider four examples 
of interest in current-day star formation: (1) a T Tauri star, (2) an embedded low-mass 
protostar, (3) an FU Orionis star, and (4) an embedded high-mass protostar. Models 1, 2, 3 
have a central star of mass 0.5 M©, and differ only in accreting at rates equal, respectively, 
to 1 X 10^^ Mq ji-^ (GuUbring et al. 1998), 2 x 10"^ Mq yr^^ (Young et al. 2003, Young & 
Evans 2005), and 2 x 10"'^ Mq yr~^ (see Popham et al. 1996, who, however, have a different 
explanation for sub-Keplerian rotation near the star than this paper). Model 4, the high- 
mass protostar, has a mass-accretion rate that is scaled relative to Model 2, the low-mass 
protostar, by their masses (see, e.g., Stauber et al. 2007), i.e., both and are taken 
to be a factor of 50 larger. In each YSO disk, we assume a standard flaring law (see, e.g., 
the dashed curve in Fig. lb of D'Alessio et al. 1999), 

A{w) = 0.1(tu/100 AU)^/^ (3-1) 



We assign tago = 3 x 10^ yr (Haisch, Lada, & Lada 2001), 1 x 10^ yr (Jijina, Myers, & 
Adams 1999), 100 yr (Herbig 1977), and 1 x 10^ yr (Osorio, Lizano, & DAUesio 1999) as 
the fiducial ages, respectively, of the T Tauri star, low-mass protostar, FU Orionis outburst, 
and the high-mass protostar. We now compute a viscous-accretion radius Ri, such that 
Md{Ru) / = tage- To ensure the approximate validity of the steady-state assumption, we 
then set R,^ = where R^ is defined as before to equal the radius that contains all the 
flux, i.e., A*(-R$) = Aq. Since Md{Ru) = M^^R^) = Manage in this formalism, the departure 
from Keplerian rotation can be computed from equation fl2-24p to equal 

1 - = (3-2) 

AgMfitage 

For protostars still building up their mass, we expect M^tagc to be comparable to M^,. On 
the other hand, for T Tauri stars or FU Orionis objects, we have M^tage small compared 
to M^. Thus, for Aq of order 4, we anticipate the departures from Keplerian rotation to be 
more substantial for T Tauri and FU Orionis stars than for low- or high-mass protostars. 

The surface density S must drop faster with than any negative power law in the 
outer parts of the disk in order to vanish, by definition, at some true outer disk edge R^- 
Therefore, unlike the assumption being made at R^,, the term represented by the right-hand 
side of equation (1-3) must change sign in the outermost parts of the disk, leading to a viscous 
movement outward of -Rd with time, carrying to large distances much of the system's angular 
momentum. Hence, the enclosed angular momentum, calculated from equation (12-261) with 
w = _R<j,, may not yield a representative estimate of the system's true total store of angular 
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momentum because it is the least centrally concentrated of the trio: magnetic flux, mass, 
and angular momentum. 

The numerical values of the relevant input parameters are now tabulated in columns 2 
through 5 of Table 2; the output bulk parameters from the theory are tabulated in columns 6 
through 11. In accordance with the discussion of §2.3 we have chosen Aq = 4 for all four cases. 
This choice results in / = 0.957 for the low-mass and high-mass protostar models, making 
magnetocentrifugally-driven, cold, disk-winds, unassisted by photo-evaporation, difficult but 
not impossible at radii much less than 100 AU for these systems, according to the criterion 
f ll-3ip for appreciable disk winds. If we had chosen a smaller value, say 2, for Aq, then / 
would have equaled 0.812, and a powerful disk wind in protostars at any radii other than 
the inner disk-edge, where magnetospheric interactions dominate, would have been almost 
as unlikely as the / = 0.658 and / = 0.386 cases tabulated for the T Tauri stars and FU 
Orionis systems with Aq = 4. 

3.1. Discussion of Bulk Properties 

In Table 2 we have chosen the numerical value of D to make = Ry (x D^^^ reasonable. 
It is informative in this regard that D = 1 works for the two protostar and FU Orionis models 
(perhaps D = 0.3 would be better), but only a relatively small value, D = 10~^'^, does as 
well for the T Tauri model (see, e.g., Andrews & Williams 2007). The enclosed disk mass 
Md{R<i>) is independent of the numerical value of D, whereas the enclosed angular momentum 
Jd{R^) varies as D"^/^. 

In the case of the FU Orionis model, we are contemplating a transient accretion event 
that has occurred during the past 100 yr, that has swept most of the mass and the magnetic 
flux to within a radius ~ -R<i> = 16.2 AU, inside of which the system is in quasi-steady 
outburst. In this case, Ry = R^ is likely to be considerably smaller that R^, where most of 
the system angular momentum may still reside. In the other cases, we think of R<p = Ry 
as being the effective "outer radius" of the system, large enough to contain the magnetic 
flux that was dragged into the system by the star formation process, but small enough so 
that the available viscosity is able to establish a quasi-steady state if the angular momentum 
contained in the system is comparable to that tabulated as Jd{R<i>)- 

Near the disk edge R^ where S becomes vanishingly small and Bz matches onto inter- 
stellar values, radial magnetic-buoyancy effects with interchange and/or Parker instabilities 
may lead to a net loss of flux from the disk. Continuing infall that brings in additional 
mass and flux from the cloud-core surroundings may counteract such tendencies. Moreover, 



Table 2. Parameters of Four Model Systems 



Object Md W D f ^ Md{R<b) Jd{R^) A(i?$) Q{R<^) 

(Mo) (Mo/yr) (yr) (AU) (Mq) (Mq AU km/s) 



T Tauri star 


0.5 


1 


X 


10-8 


3 


X 


10*^ 


10-2.5 


0.658 


298 


0.0300 


5.12 


0.480 


4.47 


Low-mass Protostar 


0.5 


2 


X 


io-<^ 


1 


X 


10^ 


1 


0.957 


318 


0.200 


51.4 


3.20 


0.381 


FU Ori 


0.5 


2 


X 


10-^ 


1 


X 


10^ 


1 


0.386 


16.5 


0.0200 


0.473 


0.320 


3.36 


High-mass Protostar 


25 


1 


X 


10-^ 


1 


X 


10^ 


1 


0.957 


1,520 


10.0 


39,700 


3.20 


0.463 



I 

I 
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equations (12-201) and (12-211) show that the outer parts of disks contain relatively more mass 
than flux, so stripping of the outer parts by edge effects, or by tidal encounters with other 
stars, or by photo-evaporation from the far ultraviolet produced by the most massive cluster 
members in dense clusters (Adams et al. 2006), cannot do much to alleviate the problem of 
the growing magnetization of the entire disk. In what follows, we ignore the complications 
that may arise from all such environmental effects. 

One could try to justify very small values of D = 10"^ to 10~^ in T Tauri models on 
the basis that only the top and bottom 10 g cm~^ of a disk (from ionization by scattered 
X-rays; see Igea & Glassgold 1999), containing a total S of 2000 or 200 g cm~^ (that apply 
roughly at 1 AU and 5 AU of the "minimum solar nebula"), are active in the accretion 
process (Gammie 1996). Such small values for the effective D are not out of the question 
if T Tauri disks have substantial "dead zones" where the ionization is too low to couple to 
magnetic fields except, possibly, for thin surface layers. The surface density ratio of active 
zone to dead zone need not be as small as 10~^'^ (see §3.2) to result in such a value for D, 
if accretion in a thin surface layer has an intrinsically smaller efficiency than fully turbulent 
MRI (see the discussion in §4.2). 

A related problem arises when we apply equation (12-101) to our T Tauri model, which 
results in the expression, a ~ 2.64(ro/AU)"'^/'^ km s~^. An isothermal sound speed of 2.64 
km s~^ corresponds to a temperature in molecular gas of ~ 1, 900 K, a value that is unlikely 
to hold even in the mid-plane at 1 AU. The difficulty arises because we took the first term 
(X A in equation (12-101) to be dominant, which requires vigorous inward accretion to sustain 
the assumed B^/ = h = 1.742 that accounts for the substantial departure from Keplerian 
rotation / = 0.658 computed in Table 2, yet we took the diffusion constant D to equal a 
paltry 10^^'^. The inconsistency disappears if we assume that the inner disks of T Tauri 
stars are dead to the MRI except for their superficial layers. With no magnetic coupling 
in the deeper layers, / would be much closer to unity in the midplane than indicated in 
Table 2, and the second term oc in equation (12-101) would dominate. We then easily 
compute that midplane temperatures at 1 AU would be closer to 600 K, and dropping as 
(cij/AU)"^/^, more in line with conventional estimates. However, we end with a non-standard 
picture where the central dead layers of T Tauri disks rotate at near-Keplerian speeds, while 
the active superficial layers tend to be very sub-Keplerian. This picture raises speculative 
conjectures that we defer to §4.2. 

The enclosed mass Md(i?$) for the T Tauri and FU Orionis models, which are inde- 
pendent of the choice of effective D (as long as it is a constant as a function of w) are 
similar to standard estimates. The disk mass of the low-mass protostar is comparable to 
those found by Jorgensen et al. (2007 and references therein) and exceed a value equal to the 
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"maximum solar nebula" when self-gravitational instabilities would need to be considered 
(Shu, Tremaine, Adams, & Ruden 1990). The disk masses of high-mass protostars are not 
observationally well studied because of their rarity and consequent larger distances. If the 
high disk masses for the two protostar models were real, gravitational instabilities could 
lead to disk fragmentation because A(-R$) and Q{R<s>) are, respectively, greater and less than 
unity in their outer disks. 

In contrast, the FU Orionis and T Tauri models are stable to disk fragmentation by both 
the magnetic and Toomre criteria, A < 1 and Q > 1. While these conclusions do depend 
somewhat on the specific choices made for D, the formation of either gaseous giant planets 
or brown dwarfs by gravitational instability at tens of AU or smaller can probably be ruled 
out in the model T Tauri system. Thermal cooling is not sufficiently rapid to help (Rafikov 
2007). Certainly, close-in brown-dwarf companions that could have been easily detected by 
the Doppler method seem difficult to produce in all our models, which is an after-the-fact 
explanation for the so-called "brown-dwarf desert" (Marcy & Butler 2000, Halwachs et al. 
2000). Our speculations in this regard are consistent with the intuitive notion that large- 
angular momentum cases are more prone to making binaries by gravitational fragmentation, 
whereas small angular-momentum cases are more prone to making planetary systems by the 
embryonic core-accumulation of solids (Lin 2006). 



3.2. Disk surface densities and magnetic fields 

Figure 2 shows E (in units of g cm~^) and Bz (in units of G), computed from equations 
f l2-8l) and fl2-7p . as functions of zu (in units of AU) for the four model systems. Hexagons 
have been placed as stop signs on the formal plots to indicate that the curves for radii larger 
than i?$ = should be ignored in any realistic apphcations. Consistent with Table 2, we 
have chosen D = 1 for the two protostars and FU Orionis, and D = 10~^'^ for the T Tauri 
model. 

In the low-mass protostar model, is 302 G, 1.09 G, and 8.74 mC at zu = 0.05 AU, 
3 AU, and 100 AU, respectively, for D = 1. In the high- mass protostar model, they are 
higher by a factor 50^/^ = 18.8 at each radius. The first value in a low-mass protostar, 302 
G at 0.05 AU, is compatible with a more-or-less smooth matching (after enhancement in the 
X-region) with the inferred field strength of the stellar magnetopause (see Shu et al. 1994), 
yielding yet another indication why disk truncation occurs inside such radii when the stellar 
field increases inward much more strongly than the extrapolated disk field. 

The second value, 1.09 G at 3 AU, is compatible with the paleomagnetism measured for 
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the chondrules of primitive meteorites, whose parent bodies are believed to originate in the 
asteroid belt (reviewed by Stacey 1976, Levy & Sonett 1978, and Cisowski & Hood 1991). 
In the X-wind model, whose predictions on this point have received strong support from the 
recent Stardust comet-sample mission (McKeegan 2006, Zolensky et al. 2006), the heating of 
the chondrule-like materials and refractory inclusions found in meteorites, and now comets, 
occurs close to the protosun (Shu, Shang, & Lee 1996; Shu et al. 1997). Nevertheless, when 
ferromagnetic chondrules are thrown out to the asteroid belt, they should not encounter 
fields that are much stronger than the inferred paleomagnetism result of 1 to 10 G (Stacey 
1976, Levy & Sonett 1978, and Cisowski & Hood 1991). 

The third value, 8.74 mG at 100 AU, offers an inviting target for Zeeman measurements 
if appropriate sources of maser emission can be found in protostellar disks. Even more 
promising for such studies, because maser emission in ring-like configurations has already 
been found (Hutawarakorn & Cohen 2005; Edris, Fuller, & Cohen 2007), are the large disk 
fields predicted for high-mass protostars if their mass- accretion rates scale anything like their 
mass even at radii of a few hundred AU. 

Consider now the FU Orionis model. At w — 0.05 AU, not far from the putative stellar 
surface, the vertical magnetic field is predicted to be 1.92 kC if D = 1, somewhat higher than 
the value for 5^ ~ 1 kC inferred from observations for the inner disk of FU Orionis itself 
(Donati et al. 2005). A four-times thicker disk, as may happen for the hot inner regions, 
would eliminate the discrepancy. Moreover, the same observations claim that the inner 
regions of the disk in FU Ori rotate at a speed 2 to 3 times lower than the Keplerian value, 
consistent with / = 0.386 in our model. A more detailed study of this system is warranted, 
but we caution that precise modeling would need to take into account the interaction of the 
magnetized accretion disk with the (squashed) stellar magnetosphere. If the mass accreted 
onto the star per FU Orionis event is Ma{Ru) ~ 0.02 M0 independent of D, as modeled 
in Table 2, then it takes tens of such events to accumulate the entire mass of the star, in 
accord with the astronomical statistics of such objects (Hartmann & Kenyon 1996, but see 
also Herbig et al. 2003). 

Figure 2 shows that the surface densities for the T Tauri model disk are 42.9 and 7.63 
g cm~^ at 1 and 10 AU, respectively, ii D = 10~^'^. These are smaller values for the planet- 
formation zones of terrestrial and giant planets than given by conventional minimum solar 
nebulas (Hayashi et al. 1985) because a comparable disk mass is spread over a much larger 
area (but see §4.2). Indeed, S at > 100 AU has dropped to such low values (below visual 
extinctions of 1 mag to the mid-plane) that any CO would be dissociated by the interstellar 
radiation field. This phenomenon and the natural steepening of E near the outer edges 
of accretion disks mentioned earlier may account for some of the larger power-law gradients 
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inferred from CO brightnesses reviewed by Dutrey, Guilloteau, & Ho (2007). Notice that the 
inferred magnetic field strength at 3 AU is 1.13 G, which is compatible with the chondritic 
values, and indeed has not changed much from the case of the low-mass protostar. The near- 
equality arises because the low-mass protostar and T Tauri models have about the same flux 
and the same radius, although their disk masses differ by a factor of 6.67 and their angular 
momenta by a factor of 10 (for the chosen values of D). Thus, Figure 2 shows that the 
surface densities differ by a factor of a httle over 6, but the two curves he almost on top 
of one another all the way out to about the same position for the two hexagons. This result 
makes graphic the point that the T Tauri disk rotates slower than the low-mass protostar 
disk because the former is more strongly magnetized relative to the disk mass. Indeed, one 
could heuristically imagine the low-mass protostar evolving into the T Tauri system if the 
excess mass and excess angular momentum could be put into an orbiting stellar companion 
without changing the magnetic flux distribution. 

With D equal to a strict constant, the surface densities corresponding to low, typical, and 
high power-law flaring, A{vj) oc w^^^, zu^^^, and zu'^^'^ are, respectively, S oc tu"^/^, ro"^/^, 
and consistent with the power-law range deduced from the thermal dust-emission of 

YSO disks (Andrews & Williams 2007), but shallower than the law E oc ccj"^/^ associated 
with conventional minimum solar nebulas (Hayashi et al. 1985). However, for T Tauri stars, 
once one allows the possibihty that D might be substantially smaller than unity because of 
physical considerations other than fully developed MRI turbulence (see below) , then there is 
no reason to think that D would be a spatial constant. On the other hand, we may do well 
to recall that the steeper, empirical, log-log slope is derived from the inferred distribution 
of solids, which may, as seems to be implied by the Comet Wild results (McKccgan 2006, 
Zolensky ct al. 2006), have been affected relative to the gas by the recycling of rock from the 
hot disk regions near the protosun to the rest of the solar nebula, as well as by the migration 
of planets. In any case, we would be the first to admit that our models do not allow for a 
straight-forward recovery of models that look like the "minimum solar nebula." Probably no 
viscous accretion disk can "succeed" in this regard (see Promang, Terqucm, & Balbus 2002). 

Vorobyov & Basu (2006) suggest that FU Orionis outbursts are associated with spiral 
gravitational instabilities in a protostcllar disk. We are sympathetic to the view that such 
self-gravitational disturbances can play a role in the early evolution of protostars that are 
still in the main infall stage that builds up the final system mass (Shu et al. 1990). We 
are however agnostic when it comes to the issue whether FU Orionis systems represent such 
early-stage objects or not. Accurate estimates of the outburst disk masses in FU Orionis 
systems - whether they are closer to "minimum" or "maximum" values - can prove to be 
observationally decisive in this debate. 
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4. HIGH AND LOW STATES OF ACCRETION 

In §3 we have presented the astrophysical case that there are high states of accretion 
where D is of order unity (FU Ori, low- and high-mass protostars) and low states where D is 
much less than unity (T Tauri stars). Indeed, even low-mass protostars (or, at least, the so- 
called Class I sources) may alternate between high- and low-states of mass accretion (White 
et al. 2007). In §4.1 that follows, we speculate that the MRI is fully developed in high states, 
and we consider mechanisms that may allow turbulent values of v and ?7, in a situation with 
dynamically strong mean fields, to achieve the saturated ratio of equation (11-241) in steady 
state. Likewise, in §4.2 that follows, we discuss the weak-turbulence conditions likely to 
prevail on scales smaller than a vertical scale height Zq in low states, focusing in particular 
on the form of MRI likely to be present when dead zones are bottlenecks to rapid disk 
accretion, with the activity being concentrated in thin surface layers (Gammie 1996). 

4.1. Magnetic loops and their dynamics 

To picture how rapid transport of matter and magnetic fluctuations across strong mean 
field lines that are anchored externally is possible, we adopt a mental image of field loops on 
a scale smaller than zq. This mental image can be given a physical correspondence in the 
equations of magnetohydrodynamics in axial symmetry, but the task becomes much harder 
if the current associated with the loop has structure in the 99 direction. We ignore this 
complication in the heuristic discussion based on a diagram (Fig. 3) that shows only one 
field loop, born of a single mean field line, that has complete freedom to move as if there were 
no constraints from neighboring field lines and other loops. Because we make no attempt 
to be quantitative except for a single order-of-magnitude calculation, this mental image can 
suggest possible interpretations without misrepresenting, hopefully, the complex, nonlinear 
dynamics of fully developed, 3-D MRI turbulence. 

Consider a process that bends, pinches, and twists a field line into a loop that eventually 
disconnects from its parent field line by resistive dissipation (bottom set of diagrams in Fig. 
3). The loop is then advected to the next set of field lines, to which it reconnects, relaxes, 
and gets into position to form another loop. To visualize what is happening in this figure, 
recall that magnetic field lines never end, but are directed continuously from point to point 
on a given line, except during reconnection, when oppositely directed fields can annihilate, 
leaving the remaining fragments to join up in a new field-line configuration. In a random 
field of fluid turbulence with a straight and uniform distribution of the mean field, a loop is 
as likely to get transported away from the star as toward it; i.e., the loops do a random walk, 
and the entire process is describable as a "diffusion" across mean field lines. The process has 
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directionality and becomes a diffusive flux when the mean field has a spatial curl defined by 
the mean accretion fiow, as drawn in the diagram. 

Note that the entire "bend-pinch-disconnect and twist- reconnect" process requires hehcal 
turbulence in three dimensions, with differential rotation providing the critical "twist" part 
of the process. Similar diagrams were drawn by Parker (1955) in his famous proposal for the 
mechanism of dynamo action. In the present context, the "twist" also implies an outward 
transport of angular momentum if it is accompanied by the shearing of the radial field 
to give an azimuthal component B^p. In this manner, the entire sequence of steps provides 
both an effective viscosity for angular momentum transport and an effective resistivity for 
matter to move from one set of field lines to another. In 2-D, "bend-pinch-disconnect" gives 
a loop that can transport angular momentum if the shear of differential rotation generates 
a from the local B^ (a "2.5 D" process). But the loop, without an additional vortical 
"twist" in the third direction (requiring an "eddy" motion out of the page, azimuthally in 
the drawing), has the wrong orientation to attach to the next set of mean field lines (top 
set of diagrams in Fig. 3 with the fourth "reconnect" step forbidden by the large red cross). 
Thus, the loop will be trapped between the thicket of nonzero mean-field lines, and it will 
eventually retrace steps 3-2-1 and merge back onto the original field line (or with other loops 
carrying the same sense of current), causing the matter to become re-attached more-or-less 
to the same field location except for the slight diffusion associated with dissipation by the 
microscopic coUisional resistivity. In such a situation, the fluctuations associated with MHD 
turbulence are probably better described as a random collection of Alfven waves rather than 
as a diffusing, merging, set of magnetic loops. 

The 3-D process appears in many MRI simulations when the plasma beta is large com- 
pared to unity (~ 100, see Appendix B). It remains to be seen if it persists in the presence 
of a mean fleld as dynamically strong as we advocate in this paper. In any case, the ran- 
dom walk of magnetic loops, carrying an associated current, can move through the thicket 
of mean field lines faster than individual particles or particulates can get knocked off one 
set of field lines by physical collisions to attach onto the next set of field lines, allowing for 
"viscous" and "resistive" diffusivities that are larger than conventional microscopic values. 
The magnetic dissipation process is sometimes described by hyper-resistivity, i.e. turbulent 
transport of current, not field, which was originally proposed to describe magnetic relaxation 
in plasmas (Strauss 1976, Diamond & Malkov 2003). Appendix A shows how the simplest 
mathematical model of a diffusion of (/^-current rather than a diffusion of 2;-field yields the 
same practical results as §2, but with a Prandtl ratio r]j/v 1 rather than rj/u <^ 1. 

In a pure- loop picture, the derivation of §2.1 really applies then to the loop dynamics 
of Figure 3. In other words, B^ of that derivation is really SB^ of the loop, with the 
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change in sign of 6Bt^ from the top to the bottom of the loop being irrelevant because the 
6B^ that is produced by shearing will have the correct compensating sign as already noted 
in the discussion of §2.1. We then assume 6B^ ~ i?+ and 5u5B^ ~ w{dfl/dzu)6zu6B^, 
with the corresponding Maxwell stress calculated from the quadratic correlation of 5B^5B^p 
assuming 5u ~ VlSvj as in §2.1. The estimate for u then goes through as before, with the 
large uncertainties in proportionality constants absorbed in JF and eventually D. Consistent 
with the discussion of the formulation of the mean- field MHD equations in §1.1, we then 
have a mathematical separation in which there is no mean B^ when averaged over z, but 
there are local fluctuations 5Bip whose correlations with 5B-^ do not average to zero. 

The estimate for t] ~ Av with A = zq/vj 1 might follow because all detached and 
non-detached loops can transport angular momentum but only a fraction of the detached 
ones have the right geometry and orientation to reconnect with mean downstream field 
lines, yielding an effective resistivity t] that is much smaller than the turbulent viscosity v. 
The exact relation (11-241) . which includes an extra factor of 3/2/^, then presumably arises 
because, in steady state, the rotation law is quasi-Keplerian and the surface density has a 
power-law index —2^ (eq. |l-22j ). The fact that the scale of the turbulent mixing length 5vj 
was left unspecified in §2.1 (in actuality, a spectrum of such scales and shapes) may give the 
problem the necessary degree of freedom to make matters come out exactly right. At some 
basic level, the macrophysics of fully developed MRI makes angular-momentum transport 
the driving energy-release mechanism behind the inward fluid drift in disk accretion. The 
formation rate, merger rate, and geometry of magnetic loops may be regulated to yield 
a turbulent resistivity, 7] = (3^4/2 J^)//, that is well below the naive Prandtl ratio rj ^ u 
because no energy source exists to cause gas to diffuse across flux tubes faster than the 
saturated value. Conversely, if t] were to fall below the level (3^4/2/^)//, the resultant pile- 
up of mean field lines waiting to diffuse inward (see LPP, whose solution for the induction 
equation remains valid independent of any implicit or explicit assumptions about fl) would 
presumably cause a shift in the numbers and kinds of loops generated until rj approaches the 
saturated level. However, more rigorous theoretical studies and/or numerical simulations are 
needed if we are to gain confidence that MRI dynamics under the circumstances envisaged in 
this paper can truly satisfy the diffusivity-ratio constraint implied by equation (11-241) with u 
given by equation (l2-6p . If such studies show that equation (ll-24p cannot be achieved, then a 
possible resolution for real systems is to alternate between trying to satisfy u = —?)u/2w and 
u = —{r]/zo)B^/Bz, making relaxation oscillations between the two conditions the real cause 
of low states and high states, with FU Orionis outbursts and their decay as the transition 
phenomenon. 

In our enthusiasm for magnetic loops, we should not forget that Alfven waves can also 
carry angular momentum, depositing it in the matter when they dissipate. If the MRI is 



-34- 



operating at maximum efficiency, it is easy to show that the frequency associated with Alfven 
waves of wavenumber scale Zq^ is comparable to fl. For larger wavenumbers (smaller scales), 
the Alfven wave frequency is larger than the natural eddy turn-over frequency Q, leading us 
to a picture of the generation of Alfven waves by the bending or twisting of proto-loops of 
scale zq that do not detach from their mean field lines, before these wavelike disturbances 
propagate, go into a free cascade, and dissipate from interactions of the type described by 
Goldreich & Sridhar's (1997) theory of MHD turbulence. It is unlikely that the competition 
with loop detachment and merging from such wave-transport effects could be adequately 
described by diffusion equations at a macro level. 



4.2. Two-dimensional MHD turbulence and layered accretion 

Creating loops of field and chopping them off from mean field lines by turbulent fiuid 
motions may not be possible when the coupling to magnetic fields is weak, and the MRI 
transport mechanism becomes confined to surface layers where the ionization level is still 
sufficiently high (the situation in T Tauri disks). The formation of the loops themselves 
becomes difficult because the restricted height practically available for z-motions may make 
the fluid effectively two-dimensional. In particular, lifting parcels of gas against their own 
weight in the ^-direction either to bend field lines or to twist them, added to the energy 
needed to stretch and pinch magnetic fields, may prove relatively difficult in thin surface 
layers compared to the same processes near the mid-plane where the vertical gravity vanishes. 
In other words, the MRI is an intrinsically 3-D instability, and it cannot operate efficiently 
in a 2-D magnetofluid except as artificial "channel flows" (see, e.g., Goodman & Xu 1994) 
where the following considerations still apply. 

In circumstances where the flow is conflned to 2-D, the turbulent resistivity rj is "quenched," 
becoming proportional to its microscopic coUisional value, although enhanced by a factor 
{SB/By when the fluctuations are large (Cattaneo & Vainshtein 1991; Gruzinov & Dia- 
mond 1994, 1966; Diamond, et al. 2005). No matter how intricately turbulence distorts 
magnetic fleld lines in the remaining two (horizontal) directions, electrically conducting par- 
ticles cannot get off the fleld lines about which they gyrate, unless they are knocked off by 
microscopic physical collisions. In 2-D MHD turbulence, an inverse cascade of squared mag- 
netic potential exists alongside the familiar energy cascade to smaller scales. This inverse 
cascade reflects a competition between the tendencies of velocity fluctuations to chop-up 
iso-contours of magnetic potential, thus producing smaller scales, and of magnetic loops to 
aggregate because of the attractive force between parallel lines of current, thus producing 
larger scales. Thus, fleld fines never get chopped up systematically into small loops that can 
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be reconnected much more quickly than the laminar dissipation of the mean fields. In layered 
accretion, therefore, the transport of mass and angular momentum are subordinate to the 
diffusion of field, and the effective turbulent v in the active surface layers may be constrained 
to be compatible with the microscopic collisional value of r/. Because the entire layer is not 
involved in the relevant diffusive processes, a formulation that integrates through the vertical 
thickness cannot do justice to the real problem which has a bimodal vertical stratification. 

At a minimum, we should consider instead a two-layer description and introduce surface 
densities, and S^, that describe respectively the columns through the upper and lower 
(active or live) surface layers and a middle (inactive or dead) layer, which sum to the total 
E of the current formulation. In this picture, the values cited for the surface densities and 
magnetic-field strengths in §3.2 probably refer more to the active layer than they do to the 
total. Although the layer thickness expressed in terms of may be more-or-less fixed by 
the (external) sources of ionization, the magnetic field strength is potentially adjustable 
as a function of w to give a constant accretion rate Md in steady state, which yields an 
advantage of such a description of layered accretion over that given originally by Gammie 
(1996). 

Our current calculations yield no constraint on the possible surface densiity of the 
inactive layer S^. Magnetic fields have freedom to move with respect to the nearly un- 
ionized gas in that rotates at near Keplerian speeds. Shear instabilities arising from 
upper and lower surfaces that rotate slower than the midplane layers could lead to breaking 
radial buoyancy waves that provide torques to redistribute the angular momentum of the 
"dead gas" (Vishniac & Diamond 1989). The coupling provided by the excitation of waves in 
dead zones has been explicitly demonstrated in the simulations by Fleming & Stone (2003) 
and Wiinsch et al. (2006). In a more simplified approach that would not attempt to resolve 
the internal structure of the upper and lower layers, the layer with surface density S^, would 
be described by the equations given in this paper, except again for a frictional term coupling 
them to Em. With included as a frictional load, the net effect would be as a variable D 
in the single-layer description of the total surface density S. In other words, v is slaved to 
T] in in the active, but geometrically thin, surface layers, and r/ is given by its collisional 
value. Then D is simply a defined quantity in the current formula 02-61) for the relationship 
between v, B^, zq, fl, and total E. Such a two-layer model, with enough microphysics to 
specify the colhsional value of in a complex, dusty, plasma, would allow us to calculate the 
variation of the effective D with w in our current one-layer description. 

Figure 4 gives the estimate by Sano et al. (2000) of the microscopic collisional value of 
rj in the mid-plane of a Hayashi model solar nebula with dust of unagglomerated interstellar 
size and ionized by Galactic cosmic rays. The collisional resistivity in the inner disk (t<7 < 3 
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AU) is consistent with the magnitude ~ 10^° cm^ s~^ cited by Shu et al. (2006) as needed for 
dynamic disk formation. The coUisional resistivity beyond ~ 10 — 20 AU is lower than the 
values needed for the low-mass protostar model (long dashes with = 1) or for the T Tauri 
model (short dashes with D = 10~^ and D = lO"'^). In the interpretation of this section, 
then, once disk formation has occurred, its accretion resistivity would need to arise from 
MRI turbulence at large radii, whereas interior to 10 — 20 AU, dead zones may be present 
and coUisional resistivities may be adequate for the needed field diffusion even in the thin 
surface layers where viscous accretion is active. 

An interesting question then arises as to what can initiate a transition between a low- 
state and a high-state of accretion. It is natural to expect the transition to originate at a 
boundary between dead zones and live zones. By definition, such boundaries are not thin 
layers in any description where vertical stratification matters. The electric fields experienced 
by charged particles forced by collisions to rotate in the dead zone relative to the magnetic 
field, at speeds characterizing the large slip between the magnetically coupled active layer 
Eg and the magnetically decoupled layer E^, may generate suprathermal particles. These 
suprathermal particles might produce the ionization that converts E^ into a better conduct- 
ing medium. Unfortunately, the existing numerical simulations of dead zones do not help us 
much in the latter regard because the acceleration of suprathermal particles requires a ki- 
netic treatment, not just a (magneto)hydrodynamic one. Moreover, the large slip is missing 
in the local simulations of Fleming & Stone (2003); and the magnetic field is missing in the 
global simulations of Wiinsch et al. (2006). Heating by the resulting enhanced accretion may 
further enhance the development of the three-dimensional, turbulence of the type with which 
we started the discussion of this section. The boundary would then eat its way radially into 
the zones that were previously dead. An interesting theoretical goal would be to see how 
this transition between low states and high states works in detail and whether an FU Orionis 
outburst begins inside-out or outsidc-in since sufficiently-ionized regions from conventional 
sources exist on both sides of normal dead zones. 



5. SUMMARY AND CONCLUSIONS 

The discussion of §4 represents our attempt to resolve the confiicts imposed by the 
following separate issues: 

a) the existence of a definite relationship between u and rj in steady state, 

b) the evidence that the common diffusion coefficient D has a value of order unity in some 
systems and much less than unity in others. 
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c) the fact that MHD turbulence has a very different character in 2-D compared to 3-D, 

d) the suggestion that r] may be hmited to have essentially its microscopic coUisional value 
in layered accretion, 

e) the difficulty of weak magnetic coupling when "dead zones" arise in YSO disks. 

Our suggestions in §4 are therefore as much a roadmap of needed future research as they are 
a catalogue of the mysteries of the present and the past. 

In retrospect, the biggest mystery concerns the most observationally well-studied disks 
associated with star formation, those in T Tauri systems. We may phrase the conundrum 
as follows. Diffusive processes cannot remove angular momentum or magnetic flux from the 
system; they can only redistribute them within the system. In a closed system, T Tauri 
stars represent an end game for viscous resistive disks whose mass steadily drains into the 
central star, but whose magnetic flux and angular momentum, inherited by gravitational 
collapse from the interstellar medium, remain more-or-less trapped in the surrounding disk. 
Such a situation must result eventually in a magnetically dominated disk. To prevent the 
residual disk from spreading to very large radii in a fixed amount of time, then demands 
inefficient diffusion (small D) . The astronomical challenge therefore becomes to explain why 
there are two physical states of accretion, an active state (protostars. FU Orionis outbursts), 
characterized by a D with a more "natural" value ~ 1, and an inactive state (T Tauri stars), 
characterized by a D <^ 1. 

The conventional assessment is that "dead zones" provide the resolution for why the 
mass- accretion rate is so low in T Tauri disks observationally, or why the effective D 'm 
our language. But if this is the correct answer, then why should D ever be as large as unity in 
other contexts, the most obvious being FU Orionis outbursts? These disks have even higher 
column densities of disk matter that can shield external sources of ionization, principally. 
X-rays and cosmic-rays. Why aren't they even more full of dead zones? We have proposed 
exploring the possibility that the high states of disk accretion correspond to the removal of 
such barriers. Perhaps a fraction of the energy released in the resistive dissipation of stressed 
fields accelerates suprathermal particles and thus provide a level of in situ ionization much 
in excess of the sources considered in conventional solar nebula models. Thus, the MRI 
mechanism, properly generalized to include the dissipation of currents generated by the 
stressing of mean fields from viscous accretion, may contain its own solution to the challenge 
posed by low ionization (sec also, Fromang et al. 2002). Heating from enhanced accretion 
may also help with the ionization of trace species such as lithium and potassium. Indeed, 
since any bootstrap mechanism allows the potential of a runaway - more ionization — > 
more coupling — > more ionization, etc. - this proposal also offers an opportunity, perhaps. 
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to understand why disk accretion in YSOs can alternate between low and high states of 
accretion, characterized generically by T Tauri stars and FU Orionis outbursts. 

The complementary difficulties noted above are related to a serious problem noticed 
by King, Pringle, & Livio (2007). The effective ccgs in a Shakura-Sunyaev prescription for 
disk viscosity has to be of order 0.1 to 0.4 to explain the empirical facts known about thin, 
fully-ionized, accretion disks in many astronomical contexts, yet the equivalent alpha from 
almost all MRI simulations to date is lower typically by one or more orders of magnitude 
(e.g., as summarized by Gammie & Johnson 2005 and modeled by Fromang et al. 2002). By 
coincidence, low values for ctss ~ 10^^ are empirically acceptable for modeling T Tauri disks 
(see, e.g., Hartmann et al. 1998), because such disks have extensive dead zones, in which 
only a small fraction of the total surface density is actively accreting. But if the estimates 
from MRI simulations were applied only to the active layers, the net effective value of ctss 
(in the sense of a combined two-layer model) would have been more hke 10~^ than 10~^. 

Among the possible resolutions mentioned by King et al. (2007) are stronger magnetic 
fields and global rather than local simulations for the MRI that occurs in realistic systems. 
The results of the present paper (see especially the discussion of §§2.1 and 4.1) strongly 
support such a resolution of the existing paradox, at least for the field of star formation. 
Most MRI simulations ignore the presence of a nonzero magnetic flux that threads through 
the disk, carried in by the process of gravitational collapse. As demonstrated in this paper, 
the presence of an externally supplied magnetic fleld makes the self-consistent dynamics 
considerably more subtle than the simplest notion of the MRI extant in the literature. In 
particular, the accretion flow generates, on either side of the mid-plane, a mean radial fleld 

from the mean vertical field because of the inward drift and the imperfect tendency 
toward field freezing. This mean radial field, whose surface value is denoted by the symbol 
B^ in this paper and whose properties can be deduced only by a global calculation that takes 
proper consideration of the vacuum conditions above (and below) the plane of the disk, sets 
the scale for turbulent fiuctuations (if MRI arises) and has three important consequences. 

First, as emphasized in §4.1, the resultant poloidal- field configuration can spawn mag- 
netic loops, possessing a radial component of the magnetic field SB^ whose amplitude is 

proportional to B^. The loop can be stretched in the azimuthal direction by the differential 
rotation in the disk to produce an azimuthal fleld SB^ that has a systematic orientation with 
respect to SB^. The correlation of 6B^ and SB^f, then exerts a Maxwell stress much larger 
than the corresponding values obtained in simulations of MRI where there is no external 
fleld i?2 to set a scale for B^. The Maxwell stress leads to angular momentum transport 
that yields the original accretion responsible for the generation of mean B^ from mean Bz. 

Second, the resultant poloidal-fleld conflguration introduces current flows that can be 
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dissipated by resistive effects. An important finding of our study, extending the work of 
LPP, is that the ratio of the effective resistivity, rj, to the turbulent viscosity, z/, must have 
a well-specified value in steady state that depends on the local aspect ratio of the disk 
(vertical thickness to radius). The exploration of the implications of this result for the 
turbulent microphysics of the problem, in particular, how the microscopic dynamics of the 
current loops can automatically adjust to the requirements of the macroscopic problem, 
needs further theoretical study, best supplemented by numerical simulation. 

Third, the resultant poloidal-field configuration produces a change in the radial force 
balance, giving a deviation from the traditionally invoked Keplerian profile. This deviation 
is not easily detectable observationally because the resulting rotation law has in steady state 
the same power-law dependence with radius as a true Kepler law, but the coefficient is 
smaller. Thus, even if it were present, observers would tend to attribute the result to the 
mass of the central object being smaller than its actual value, or to the disk being inclined 
by a lesser amount than in reality. Nevertheless, it would be illuminating to find such an 
effect in the YSO disks that have the smallest mass-to-fiux ratios. Indeed, a deduction of 
sub-Keplerian rotation of the disk may already have been made in the case of FU Orionis 
(Donati et al. 2005), but the proper interpretation of the phenomenon in this case may be 
complicated by the interaction of the disk field with the imperfectly squashed stellar fields of 
the central object. Finally, it has not escaped our attention that significant departures from 
true Keplerian rotation of a YSO disk may have important consequences for the problems of 
binary-star and planetary-system formation and evolution, particularly with regard to the 
difficult issues of orbit migration and eccentricity pumping (e.g., Goldreich & Sari 2004). 
The inward drift of solids from sub-Keplerian regions into the dead zones of the problem, 
which rotate more nearly at Keplerian speeds, may give an additional reason to focus on 
the importance of dead zones for the problems of planetesimal and planet formation (e.g., 
Youdin & Shu 2002, Pudritz & Matsumura 2004). 
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support from CONACyT 48901 and PAPIIT-UNAM IN106107; AG, from NSF Grant AST- 
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A. ALTERNATIVE FORMULATIONS FOR TURBULENT DIFFUSION 

A turbulent, magnetized medium may not behave in the assumed model fashion, with 
diffusive fluxes proportional to a scalar diffusivity times the rate of spatial variation of 
the mean quantity that is being spread (Pick's law). For example, Pessah et al. (2007) 
perform a quasi-linear analysis of the MRI with third-order closure in a simple shearing-box 
geometry. They claim that the turbulent viscous stress 11^^? depends on the rate of shear 
by a non-Pickian power p that is different from 1. If we denote {w/Q)dQ/dw by —S, then 
their model yields H^^jip = —OTiQS^ with p between 3 and 4. However, in a field of quasi- 
Keplerian rotation where Q oc zu~'^/'^, S is simply the number 3/2. Thus, the Pessah et 
al. formalism satisfies, in practice, the usual Newtonian relationship, n^i^ = u'EwdQ/dvu, 
where z/ = (3/2)^~^i>. This simple transformation allows the translation of all of our results 
into the language of Pessah et al. 

A similar remark applies to the induction equation fll-41) . Instead of the diffusion of 
vertical field, some mean- field formulations of turbulent MHD (see §4.1) envisage the diffusion 
of tangential current (per unit length). Prom Ampere's law, = (c/27r)_B+, and the diffusion 
of mean current is equivalent to the diffusion of the mean radial field. Then, instead of 
equation (11-41) . we might postulate a diffusion equation of the form 



where rjj is the diffusivity for the turbulent diffusion of current. In steady-state, radial force 
balance in a field of quasi-Keplerian rotation will still require = hBz oc w^'^^^'^'^^/^ (see 
§§1.3 and 2). In this case, the above equation is equivalent to equation (11-41) if we identify 



where we have used equation fll-24p . For n = 1/4, the relationship between t]j and u is 
then rij = 1.09 z/, a result that we might call the "Prandtl hypothesis" because it differs 
only slightly from the naive guess that diffusion processes in a turbulent medium have equal 
steady-state diffusivities (cf. the discussion in LPP). Thus, what seems more relevant than 
the specific turbulent diffusivities are the turbulent fluxes, and how those fluxes relate to the 
spatial derivatives of mean-flow quantities. 




(Al) 




(A2) 



B. THE VALUE OF D FROM MRI SIMULATIONS 



In Table 3, we list the ratio of magnetic stress to magnetic energy density from Table 
1 of Stone et al. (1996; SHGB) and from Table 1 of Miller & Stone (2000; MS). Double 
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angle-brackets indicate time and volume averages (mostly over two scale heights) for the two 
smallest values of the initial mid-plane beta parameter, /3(0), in these papers. SHGB use 
periodic boundary conditions in all three directions and a vertical box height three times 
the initial pressure scale height H . MS employ an outgoing boundary condition at the top 
of a box of height bH. In the Table, two averages are included for the results of MS, which 
apply to the regions \z\ < 2H (fourth row) and \z\ > 2H (fifth row). Empty entries occur 
when the required information is not given in the original papers. 

In the simulations, and By contain only fluctuating components, whereas Bz and 
total B contain both mean and fluctuating components, with the mean components being 
systematically destroyed by numerical reconnection as the simulations proceed because the 
starting distribution of B^ alternates in sign radially. If we equate the tangential Maxwell 
stress per unit circumferential length (which dominates the corresponding expressions for the 
Reynolds stress) to an equivalent viscous stress in the usual manner, we obtain the following 
expression for u for the case of quasi-Keplerian rotation, 

- {{B^By)) (Bl) 



Introducing a factor of {{B^)) on top and bottom to express the relevant ratios in the form 
of Table 3, we can now identify the coefficient D in equation (12-61) as 

1 {{B^By)) {{B^)) 

6vr {{B^)) m))' ^ ^ 



The results in Table 3 give values of D that differ by an order of magnitude, indicating 
that the MRI simulations are sensitive to the assumed boundary conditions and to the size 
of the computational box. As a formal result, the last column of Table 3 seems to suggest 
that D is, at best, 0.2. However, as emphasized in the text, simulations with net magnetic 
fiux equal to zero do not correspond to the situation of interest for our study. 



Table 3. Estimates of D 



Run 


m 


mBy))/{{B^)) 


{{B^))/{{Bl)) 


D 


SHGB IZl 


100 


0.145 


26.7 


0.2 


SHGB IZ3 


25 


0.139 






MS ZN2 


100 


0.0716/0.0958 






MS ZNl 


25 


0.0111/0.00586 


40/26 


0.03/0.01 
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C. VERTICAL STRUCTURE OF STRONGLY MAGNETIZED DISKS 

Inside a thin disk, the condition of vertical hydrostatic equihbrium reads 

oz \ Svr / 

where P is the thermal gas pressure and is a function of z equal to at 2; = and to 
at the disk's surface, whereas B^ may be taken to be a constant over the same range of 
z. If we integrate equation (ICll) from ^ = to the surface under the boundary conditions 
that -P(O) = a'^'E/2zQ (which defines what we mean by a^) and P = at the disk's surface 
(which defines what we mean by the surface), we get 



Svr 2zq 4tz73 ' 



(C2) 



where we have defined Zq by requiring the integral oi zp from 2; = to the surface of the 
disk yield (2;o/2)(S/2). With = /^i?^, the equation of radial force balance (11-181) reads 



Mi = n _ f2^GM. 

27rE ^ ^ ' w2 



f)^. (C3) 



Elimination of Bl from equations ( lC2p and ( 1C3I) . with zq = Azu, then gives equation f l2-10p . 
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Fig. 1. — Schematic diagram of accretion flow in a disk threaded by magnetic flux accumu- 
lated by the process of star formation. 
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Fig. 2. — The vertical component of the magnetic field (solid curves) and the surface 
density S (dashed curves) plotted against the radius zu in the steady-state disks of the 
four models of Table 2. The hexagons mark the location where Ru = in the T Tauri 
(red curves), low- mass protostar (black curves), FU Orionis (blue curves), and high-mass 
protostar (green curves) models. The slopes are —11/8 and —3/4, respectively, for logB^ 
and logS versus \ogzu. 
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Fig. 3. — Schematic diagram of scenarios by which field loops are created by magnetohydro- 
dynamic turbulence when the mean field is strong: (top) in 2-D by stretch, pinch, disconnect 
and (bottom) in 3-D by stretch, pinch, disconnect and twist, reconnect, relax. The depiction 
is the meridional plane {zu,z), except for the twist indicated by the block arrow, which oc- 
curs out of the plane of the paper in the ip direction because of differential rotation. Notice 
the bias for forming the loop on the side closer to the star because of the accretion flow. 
This bias causes the diffusive flux to flow in the correct direction relative to the curl of mean 
B. Because the loop in the top diagram does not experience the twist operation, it has the 
wrong orientation to reconnect with the mean fleld downstream of the mean accretion flow 
since the flelds point up on both sides of the target contact point. The twist in the bottom 
diagram gets the flelds oriented in opposite directions at the target reconnection point, which 
results in the green "yes" sign to proceed to steps 4 and 5. 
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Fig. 4. — Comparison of coUisional resistivity (solid curve) applicable to the Sano et al. 
(2000) model of the minimum solar nebula, which includes the effects of cosmic grains of 
typical interstellar size (a = 0.1 /xm) with the required turbulent values in the models of §3 
of a low-mass protostar (long dashes) with D — 1 and a T Tauri star (short dashes) with 
D = 10^^ (upper curve) and D = 10^^ (lower curve). The hexagons mark the corresponding 
locations of i?$ where the disk holds the trapped flux corresponding to a dimensionless 
mass-to-flux ratio Aq = 4. 



